Two- Channel Kondo Model as a Fixed Point of Local 
Electron-Phonon Coupling System 

Hiroaki KUSUNOSE and Kazumasa MiYAKE 
Department of Material Physics, Faculty of Engineering Science, 
Osaka University, Toyonaka 560 

Abstract 

It is shown on the basis of the multipHcative renormaHzation-group method of 
two-loop order that the low-energy effective Hamiltonian of a strongly coupled 
local electron-phonon system is mapped to the two-channel Kondo model. A 
phonon is treated as an Einstein oscillator with restricted Hilbert space such 
that up to one-phonon process is taken into account. By eliminating the 
high energy process of conduction electrons, it is shown that a certain class of 
couplings between ion vibrations and conduction electrons is selectively grown 
up. As a result the system is reduced to the two-channel Kondo model. 
The crossover temperature Tk and the renormalized phonon frequency 
are expressed in terms of the mass ratio m/M, m and M being the mass 
of electron and ion, and the electron-phonon coupling g/D, D being half the 
bandwidth of conduction electrons. The anomalous behaviors associated with 
this renormalization can be mesuarable if the condition Tk > is fulfilled. 
It is demonstrated that such condition is satisfied when g/D is sufficiently 
large but in a realistic range. 
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I. INTRODUCTION 



It has been recognized for these decades that a certain class of A15 compounds have 
puzzhng physical properties One of the most important feature of these compounds is 
the large anharmonicity of ionic oscillations, which is exhibited, for example, as extremely 
small Debye- Waller factor at zero temperature limit and anomalously large resistivity of 
the loffe-Regel limit at room temperature, where the mean free path of conduction electron 
is comparable to the lattice constant 0. These compounds also show anomalously strong 
temperature dependence of electronic properties such as the magnetic susceptibility x 1 the 
Korringa constant 1/TiT (1/Ti being the nuclear relaxation rate) [^, and so on 0. These 
anomalies can be explained by assuming an existence of extremely sharp peak in the density 
of states, of the order of 10 ~ 10^ K, around the Fermi level. Such a sharp peak should arise 
from the many-body effect because the position of the peak always stays right at the Fermi 
level independently of either components of compounds or degrees of stoichiometry p|. 

Similar anomalies are exhibited in "heavy fermions" where the conduction electrons and 
an array of localized spins form a coherent quasiparticle band with extremely narrow width 
below the coherent temperature Tq. On the basis of this observation, the analogy between 
these two systems was pointed out In A15 compounds, the origin of the "spin" degrees 
of freedom is attributed to doubly degenerate stable positions of ion in a possible double- 
well structure of ionic potential at A site of A3B compounds 0. Similar idea was proposed 
as a model for a layered compound 2if-TaSe2 [H- Such an anharmonicity of local-phonon 
can arise through strong electron-phonon coupling [^P,^, which causes screening of ionic 
potential the electron clouds feel and in turn decreases the restoring force for the ionic 
displacement leading to softening of ion oscillations. If the coupling is strong enough, a 
simple harmonic potential for the ion may deform into a double- well [0,^. 

In order that such system is renormalized into the Kondo model as in heavy fermions, 
we have to keep the electron-ion coupling containing at least quadratic term with respect to 
the ion displacement 0. This is because it is such a term that gives the hopping between 
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the two stable positions in the double-well corresponding to the spin-flip scattering in the 
Kondo problem which is indispensable to obtain the Kondo scaling. Along such a scenario, 
a crude and qualitative description of strongly coupled electron-phonon system was given 
previously 0. 

The purpose of this paper is to investigate this problem more systematically without 
assuming the existense of double-well but on the model of an Einstein oscillator interacting 
with electron gas. By restricting the Hilbert space of phonon up to the one-phonon process, 
the Hamiltonian is transformed to the form which can be treated by the formalism of Vladar 



and Zawadowski [T^. It is shown that the Hamiltonian would be always renormalized into 
the form of Kondo model if the scaling step were not interrupted by the renormalized level 
splitting of ion oscillations which plays a role of an external pseudo-magnetic field acting 
on the pseudo-spin. In order that the "Kondo regime" is realized, the energy scale Tk 
characterizing the "Kondo renormalization" is larger than the renormalized level splitting 
A^'. Then it is shown that the condition, Tk > A^', is satisfied in the system with large but 
realistic value of electron-ion coupling. 

In the region where the Kondo model is applicable, the effective Hamiltonian can be 
regarded as that of the two-channel Kondo model [^T|-|T^ because the conduction electrons 
have the real-spin degrees of freedom other than that of the pseudo-spin which becomes 
explicit through the Vladar- Zawadowski renormalization. Here the real-spin degrees of free- 
dom play a role of the channel there. It has been well recognized that the two-channel 
Kondo model exhibits the non-Fermi liquid fixed point in contrast with the single- channel 
Kondo model. Indeed, the two-channel model has been fully solved by variety of methods 
T^-pO[. Quite recently, it has been shown that the susceptibility of real-spin, or channel. 



does exhibit the non-Fermi liquid behavior if the perturbations breaking the particle-hole 
symmetry, such as the Coulomb repulsion among conduction electrons or potential scatter- 



ing, are taken into account |2T|]. This may potentially give an explanation for the anomalous 
temperature dependence of the spin susceptibility x iii A15 compound VaSi 0]. 

This paper is oraganized as follows. In §2 we present a local phonon model interacting 



with conduction electrons and derive a simplified model describing the low-energy physics. 
Next in §3, we derive the scaling equations following the argument by Vladar and Zawad- 



owski |Ty] on the basis of the formalism of multiplicative renormalization-group [^,^. In 
§4 we solve the scaling equations and determine the crossover temperature Tk and the renor- 
malized level splitting A^. On this basis we discuss the possibility of observing the anomaly 
associated with the two-channel Kondo model. In the final section, we summarize the results 
and discuss their implications, especially the relation with the so-called Migdal approxima- 
tion which states the absence of the vertex correction for electron-phonon coupling and of 
the renormalization of the spin susceptiblility p4|,p5 . 



II. MODEL HAMILTONIAN 



A. Local Phonon Model 



We consider an Einstein oscillator interacting with conduction electrons. This simulates 
an optical phonon in A15 compounds where the transition metal ions maintain almost the 
atomic nature. In such a situation, an electrostatic restoring force tends to be screened by 
the cloud of electrons. This increases the anharmonicity of ionic oscillations and can make 
the adiabatic potential for ion displacement be even double-well like if the electron-phonon 
coupling is sufficiently strong. In order to verify such a scenario, we start with a Hamiltonian 
given as follows: 

H = Hc\ + -ffph + -f^ep, (1) 

^el = Xl^kflLaka, (2) 
k,(T 

i/ph = n (b^b + , (3) 

Hep= J2 J ^"^^(^ ~ Q^) exp[i(k' - k) ■ r]aj^^ak'a, (4) 



kk', 



where a|^^ and 6^ denote the creation operator for the conduction electron with the wave 
vector k and the spin a and for the Einstein phonon with the energy Q, respectively. C,k is 



the kinetic energy measured from the Fermi level. 

The electrostatic potential v in (H) is assumed to be contact type because its range is 
expected to be very short of the order of atomic radius as mentined above: 

v{y - Qz) -g5{Y - Qz), (5) 

where Q is the displacement of the ion along z-axis and (? is a coupling constant, which 
is positive and of the order of a bandwidth because its origin is the Coulomb attraction 
between electons and the ion. Thus the electron-phonon interaction (1) is reduced to 

-^ep = -9 5Z IZexp[?(fc; - k,)Q]al^a]^,„. (6) 

k,k' cr 



B. Simplification of the Model 

As it will be shown below, the interaction between ionic vibrations and electrons near 
the Fermi level increases logarithmically as eliminating high-energy processes. Namely, the 
electrons in the vicinity of the Fermi level is crucial; so that the polarization of conduc- 
tion electrons, which is expressed in terms of directional dependence of the wave vector k 
near the Fermi level, plays a crucial role. Therefore, we first introduce the spherical wave 
representation for the creation operator of conduction electrons as follows: 



klma \ / 

'OTT 



/ dkYUk)al^. (7) 



where Yim is the spherical harmonics and R denotes the radius of the system. With the use 
of a linearized dispersion for conduction electrons, S,k ~ kp{k — kp)/m, kp being the Fermi 
wavenumber and m the mass of conduction electron, H^i given by @ is reduced to 

Hei = J dkkali^^akima, (8) 

where the bandwidth D of conduction electrons is given as 

D = ^. (9) 

mvr 



Secondly, we make a simplification of the phonon part Hp^, (|]). The displacement Q 
and its canonical momentum P are represented in terms of the phonon operators, 6^ and b, 
as follows: 

Q = g(b + b^), (10) 

p = ^(;-*'). (11) 



where q = y l/2Mf2, M being the mass of the ion. Since the low-energy phonon states are 
important, we restrict the Hilbert space of the phonon in such a way that only the states 
with n =< b'^b >= or 1 are included in the low-energy effective (fixed point) Hamiltonian. 
Then, Q, (p!0D, and P, (pT]), are represented in this restricted Hilbert space as 

Q = qT. bir^^'bw, Q' = q'Y: bU26nn' - r„^„0&n', (12) 

nn' nn' 

P=Yi:birL'bn', P^ = ^E&U25nn'-C)&n', (13) 
^ nn' " nn' 

where r* [i = x, y, z) is the z-th component of the Pauli operator, and b\^ is the pseudo- 
fermion operator creating the n-phonon state so that |r can be regarded as the pseudo-spin 
corresponding to the phonon degrees of freedom. 

It is noted that {Qf ^ (Q^) and {Pf ^ (P^) since (Q^) and (P^) include virtual 
2-phonon process. Indeed, the matrix elements of are given as follows: 

< 0|Q2|o >=< o|Q|l >< 1|Q|0 >= (14) 

< l\Q^\l >=< l\Q\0 >< 0\Q\1 > + < 1\Q\2 >< 2\Q\1 >= ^q^, (15) 

< 1|Q2|0 >=< 0|Q2|i >= 0, (16) 

where \n > denotes the ra-phonon state. These matrix elements are equivalent to the ex- 
pression (0). One can show that the same arguments hold for the operator P^. 

For the later discussions, it is more convenient to introduce the alternative basis for the 
representation of phonon states as follows: 

b\ = ^{bl-b{), b\ = =^ibl + b{). (17) 



Then the non-vanishing matrix elements of Q are given as 

<TIQIT>=-g, <i\Q\[>=+q, (18) 

where | t> denotes the state 6||vac>, |vac> being the vacuum state, and so on. 
In this representation, Eqs. (0) and (|13|) are transformed as 

Q = -?Efe/.&/^, = q'Y:hi{26^p + Tl^)hp, (19) 

T i T J. 

P = Ya^birl^h,, = ^,Y.bi{25^p + r:,)h,, (20) 

" a/3 " a/3 

and ifph, (0), can be written in the restricted Hilbert space as 

T i 

i7ph = ^fiEfeV/3 + ^- (21) 

a/3 

Thirdly, we simplify the electron-phonon interaction Hep., (^. Since the exponent of 
exp[i(k' — k) ■ r] in (Q) is of the order of k-pQ ~ k-pq ~ [m/MY^^ <^ 1, we expand the 
exponential with respect to kpQ up to second order. Then, using the expressions (P^, we 
obtain 



afS 



+ 0{{kpQf), (22) 



Hep = 4a«k'a 
kk' o- al3 

V'^^, = -\{K-kfklq\ \f.=0, V^^,, = - 4)M, (23) 
\/A = 1 _ (4 _ KfkW- (24) 

It is noted that the interaction V^^,, (p3D, arises from Q^-term and gives the pseudo-spin flip 



scattering which is the heart of the Kondo effect . 



In order to treat the problem in the spherical representation, we introduce the interaction 
matrices VJJ, defined by 

Vl^ = 47r 5: z'''%,{k)Y*,{k')Vi,. (25) 

It is noted that the phonon vibration along the 2r-axis interacts only with the component of 
m = 0. It is convenient to introduce the basis of conduction electron corresponding to those 
of the phonon, (p!?!), as follows: 



-j- 


! 1 

"^{^kOOa + %10o-)' 


(26) 


t _ 
kla — 


"^i^kOOa ~ '^IlOo-)' 


(27) 


t _ 
kdcr 




(28) 



With the use of (0) and (pSj), the simphfied Hamiltonian is given as 

H/D = Hc\ + i/ph + Hep, 
^ei = I] I] / dkkali^akia, 

Hep= dk dk'J2J2J2Yl o.kiAi'(^k'i'abiT'^^b 



(29) 
(30) 

(31) 
(32) 



where = —gVi]/ are the electron-phonon couphngs which are non-dimensional and A* is 
a fictitious magnetic field acting on the pseudo-spin ^r. The explicit forms of vy and A* in 



the new basis are given as follows: 



-X ^ 9 1 2 2 
V = KpO 

3D 



^ 1 -1/^/TO 

1 -i/VTo 
y -i/Vio -i/VTo 



(33) 



^ o^ 



1 n 

2D' 



/\y = A' 



V 

0. 



0-10 




(34) 



(35) 



It is noted that the pseudo-spin r appears in the Hamiltonian ( p9D representing the phonon 
degree of freedom, while it does not in the electronic degrees of freedom at the biginning 
but it is induced by the renormalization group evolution as discussed below. It is also noted 
that the potential scattering is neglected in (p9D, although it has been recently suggested 
that the potential scattering is relevant to discuss the magnetic property on the basis of the 
pseudo-spin model . 
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III. SCALING EQUATIONS 



The Hamiltonian (^9|) involves the logarithmic divergence at the Fermi level due to the 
screening of local internal degrees of freedom by the conduction electons as discussed by 



Vladar and Zawadowski [10|. Since we are intersted in the physics near the Fermi level, we 
eliminate the high energy process in the sense of the renormalization group evolution. To 
this end, we apply the multiplicative renormalization-group P^IES formalism with help of 



the Abrikosov pseudo-fermion representation of the pseudo-spin 



A. Multiplicative Renormalization-Group Transformation 

First, we discuss the formalism for the multiplicative Renormalization-group transfor- 
mation on the basis of the perturbation theory. 

It is useful to introduce the Matsubara Green functions for the electron and the pseudo- 
fermion defined as 

G = (36) 

and 

= -■ T ^ Aj i ; (37) 

respectively, where u and e denote the Matsubara frequency. The electron self-energy con- 
tains a closed pseudo-fermion loop and it tends to zero as A — >■ oo; thus, the electron self- 
energy is ignored. It is noted that Q and the self-energy S for pseudo-fermion are matrices 
in the pseudo-spin space. 

It is assumed that the system with a reduced bandwidth D' behaves at low energies in 
a way similar to that of original system with the bandwidth D, if the couplings vy and the 
phonon energy A* are modified apropriately. The difference between the original and the 
scaled Green function and the vertex function can be given by multiplicative factors Z2 and 
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Zli,, respectively. The multiplicative factor for the conduction electron is given as Zi = 1. 
The multiplicative renormalization-group transformation can be given by 

Gap 4, A'^) = Z, vl^ G^p v\,, A^) , (38) 



llv 



D' 



UJ 



Ya'{^,vl,), (39) 



and 



4 = Z-'Zi,vl„ (40) 
where Yu, is the normalized vertex related to the vertex T^, as 

^i'=7wvii'- (41) 

The new scaled couplings v^, and parameters A'* are labeled by prime. An alternative 
multiplicative renormalization relation for the vertex is given by multiplying Eq. (^9|) by v'll, 
and by inserting the Eq. (^) and it is obtained as 

(^,,^u) = Z^-' (^,vu] (^,vi^ . (42) 



The multiplicative factor for A is not simple and it will be given only along the detailed 
calculation. It is important to note that the multiplicative factors depend only on the relative 
change of the bandwidth D'/D and on the couphngs vy. 



In order to construct the renormalization-group transformations given by Eqs. (38) - 
(^), the perturbation theory is applied. Therefore the result to be derived is not valid in 
the region of strong coupling. According to the perturbation theory, the scaling equations 
is given as 

^vl,=(3{vi,), (43) 

and 

^A^ = /(^;,V,A^), (44) 

where x = ln{D/D'). The irrelevant part of the phase space for the conduction electrons 
will be eliminated by integrating Eqs. (^) and (0), but the procedure can be applied only 
as far as x < min[ln(D/fcB^); ln(l/| A'|)]. 
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B. Derivation of Scaling Equations 



In order to construct the scaling equations, it is convenient to use a new basis giving 
7^ 0, which is obtained by the rotation by the angle 7r/2 around the y-axis, instead of the 

basis where 7^ and A^ = 0. In this new rotated representation, the couplings and the 

fictitious magnetic field are related as = , = v'^ , = —v^, and A^ = A^, A^ = A^, 

A^ = — A^', respectively. 

The vertex corrections of first and second order are shown in Fig. 1(a), 1(b) and Fig. 2, 

respectively. Their analytical expressions are given as follows: 



jk 



D 



\UJ\ 



(45) 



and 



rf)=n$:[2Tr(t;V>/,-Tr( 



\uj\ 



(46) 



where n represents a degree of freedom of spins of the conduction electron, i.e. n = 2. 
The self-energy shown in Fig. 3 is given by 

D 



UJ 



Tt{vV)6^p - {2Tr(t;V)A^ - Tr{v'v')A^}f- 



In 



\UJ\ 



(47) 



It is noted that the self-energy contains off-diagonal terms which are proportional to r^^ and 



r^^. Let us define A' such that the renormalization factor Z2 is independent of A*. Explicit 



expresssions of A'* and Z2 are given by 



A'' = A' + 2nY [Tt{v'v^)A^ - Tt{v^v^)A^' 



In 



D 



D 



Z2 = l + n^Tr(i;V)ln— . 



(48) 
(49) 



It will be shown later that A^ and A^ can always be eliminated by a rotation around the 
y-axis and x-axis in the pseudo-spin space. In the case A* = A^^j^, Eqs. (|^) and (|i8| ) are 
rewritten as 
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-n 



to 



X 



D 



In 



D 



\UJ\ 



A' = 2nA'TT{vV) In — , {i = x,y), 
A'" = A" 1 - 2nTr(t;^i;^ + In ^ 



(50) 

(51) 
(52) 



respectively. 

The renormalized coupling v'li, is obtained by Eqs. (^2|) , (^51) , (^61) and (^) as follows: 



= ^;,V - 2* E(^'^')"'e*'' In ^ - 2n E [Tr(t;^t;^>,V - Mv'v'H] In ^. (53) 



Although the renormalization-group transformation generates the parameters A'^ and 
A'^, they can be eliminated by the rotation around the x- and y-axes. The angles of the 
rotations are 



ax = 2nTr(t>^f ^) In — , 
ay = 2nTr(i;^r)ln— . 



Then, A' is transformed to 



(54) 
(55) 



A""" = A"^ = 0, 



A" = A^ 



D 



1 - 2nTr(w^w^ + vH^) In — 



Furthermore, these rotations modify the couphngs v'li, to v"li,, 



^ w 



D 



v'w - '^nYY{v'v')vly In — , (z = x, y) 



D 



v"l, = v'l, + 2n [1:t{v^v^)v^, + Tr(t;^r)t;,^,] In - 



(56) 
(57) 



(58) 
(59) 



Combining the relations (p3D , (|58|) , and (0), we obtain the scaling equations for the 
coupling w's as 

d 



dx 



jk 



E{Tr(t;^t;^>,V-Tr(t;V>/,} 



12 



+ TT(v'v')vf„ 



d 

dx 



jk 



(60) 



(61) 



and the scaling equation for the fictitious magnetic field as 



dx 



(62) 



These scaling equations can be represented in the original basis of pseudo-spin, i.e., in 
terms of f 's and A's without bar, if the pseudo-spin axis is rotated back around the y-axis 
by the angle —tx/2. Then, the scaling equations are given as follows: 

d 



dx 



d 



jk j 



Ti{v''v')vfi, - Tr(t;"t;^)4 



jk 



{I = y,z) 



d 

dx 



(63) 

(64) 
(65) 



IV. RENORMALIZATION-GROUP EVOLUTIONS 

In this section, solving Eqs. (|53D-(p3|), we determine the characteristic temperature Tk, 
which characterizes the crossover between weak and strong coupling regime, and the renor- 
malized first excited energy of phonon A^'(x), below which the the renormalization-group 
transformation cannot be proceeded further. Then, the typical cases of interest will be 
discussed. 
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A. Solution of Scaling Equations 



First, we discuss the case \vfi,\, \vfi,\ <^ v^i,, which appears at the initial stage of 
renormahzation-group evolution of the present problem. In this case, we can linearize the 
scaling equations ( |63D and (34) with respect to vfi, and vfi,. Then, the linearized version of 
the scaling equations, ( |B5D and (p^), are given as follows: 

-^vfi, = -2i[v^, v']u> - 2nTT{v'v')vfi, + 4nTr(t;"t;")t;f;„ (66) 
-^vfi, = -2i[v\ v%v - 2nTT{v'v')vy, + 2nTr(t;V)t;f;,, (67) 

These equations are valid in the region of x where vfi,{x), vfi,{x) <^ ^f//(0). 

We have verified by numerical calculations that the last term in Eqs. ( pBD and ( pTj ) are 
negligible PB[, although we do not present it here. Then, if we choose the representation 



where v^, is diagonal, i.e. vfi, = vfSw, the scaling equations above can be solved separately 
with the boundary condition ff;,(0) = in the following forms 

vu'ix) = vUO) cosh(2K(0) - vf{0)]x) exp{-2nTT[v\0)v\0)]x), (69) 
vy,{x) = ivfi,{0) sinh(2[wf,(0) - vf{0)]x) exp(-2nTr[t;^(0)t;^(0)]x), (70) 
vfAx)=vfmi'. (71) 

This result indicates that the couplings vfi, and v^, with the combination of / and /', 
for which \vf,{0) — (0)| takes the largest value, increases most rapidly as x increases by 
the renormalizaton-group step. Therefore, we are left with only the 2x2 subspace in the 
3x3 space of matrices fz/''s; i.e, we can effectively describe the polarization degrees of 
freeedom of conduction electrons also by the pseudo-spin of 1/2, i.e., / =1 and |, after 
the renormalization evolution is proceeded enough. In the subspace, the solutions can be 
rewritten in terms of the Pauli matix ay for conduction electron as vli,{x) = f*(x)cr;;,. Using 
these solutions, we can confirm that the last term in Eqs. (|66D and (^) can be neglected. 
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It is noed that the renormahzation as above arises only if we take into account the 
electron-phonon (ion) couphng V^^,, (p3D, including at least quadratic term with respect 
to the ion displacement because such a term gives the pseudo-spin-flip scattering with 
conduction electrons which is the heart of Kondo renormahzation. 

The 2x2 subspace is well-defined in a region x > Xi, where Xi is determined by the 
conditions 



t;^(xi) =vy{xi) = -v^{0)exp{4v'{0)[l-nv'{0)]xi), 



(72) 
(73) 



In the region x > Xi, the pseudo-spin degrees of freedom for conduction electrons becomes 
explicit as the anharmonicity of ion vibrations grows. Thus our electron-phonon model can 
be mapped to an anisotropic two-channel Kondo model. 

Next, we discuss the region v^{x) = f^(x) ~ v^{x). In this region, our scaling equations 
can be simplified as 



^v' = 4:{v^Y[l - 2nv% 
ax 

^A" = -4nA"[(i;")2 + (t;")=']. 



(74) 
(75) 
(76) 



As will be shown in Appendix, the solutions of these scaling equations, (|71D-([7BD, are given 
by 



77 1 

x = ln[v%0)v'm + — In 



4f^(0) 
f^(0) 



+ 



^f^(O) 



In 



v'^ix) +f^(0) 



(77) 



which is an implicit equation for v^{x) and v^{x) is also given in terms of v^{x) by the 
following relation 



V [XI 



[v^ix)f 



. 1 — 2nv^(x) 



1/2 



(78) 



1 - 2nv'{0) J 

The crossover temperature Tk is given so as the condition v^{x)/v^{0) ^ 1 is satisfied for 
D' < ksT^, where the last term in Eq. ( [77D can be ignored. Then solving Eq. ( [77| ) without 
the last term, we obtain Tk as follows: 
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^■<-^ = (l7i))""'°*'"^(«>"'(°M"'* (™) 

In the crossover region, the second order terms of the /5-function, i.e., the right-hand side of 
the scahng equations, become comparable to that of first order. Thus the higher order terms 
of renormahzation must be taken into account for lower energy process in gerneral. Namely, 
the perturbational renormalization-group theory can not be applied below Tk. However, it 
has been known that the two-channel Kondo model has the non-trivial fixed-point which is 
located in the intermediate coupling regime given as 

2n ^ ' 

So, does not diverge even below Tk in contrast with the conventional single channel 
Kondo problem. In this sense, Tk gives the energy scale where the crossover from weak 
to intermediate coupling occurs. Around the fixed point, (|80D, the effectve Hamiltonian 
describing the low-energy physics is expected to be given by that of an isotropic two-channel 
Kondo model. 



B. Renormalization of Phonon Excitation 

The first-excited energy of ion vibrations, A^, is renormalized by the scaling equation 
(|65|): A"^ softens as the renormalization process is proceeded because Tr(f^f^ + v^v^) > 0. 
This arises from the self-energy ( pOf ) and represents the physical process of screening 
of spring constant of ion vibration by the conduction electrons. Such an effect causes 
the anharmonicity of the ion vibrations logarithmically so long as the condition x < 
mm[ln{D/kBT),ln{l/A^{x))] is not broken, even if the bare couplings are weak. 

The softening of A^' is caused mainly by the coupling with v^v^ in Eq. (p3|), because 
is the largest coupling through out the renormalization steps. So, we estimate it using the 
scaling equation ( [76D without the term (f^')^, i.e., 

^ ln(A") = -4n(t;^)l (81) 
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As will be shown in Appendix, the solution of this equation is given by 

A^(x) = A^(0)[1 - 2nv^{x)Y/' [2^)) ' ^^^^ 

Since the energy A^(x) is required to exchange, or flip, the pseudo-spin dynamically, it 
looks as if the pseudo-spin associated with phonon degrees of freedom is quenched for the 
scattering with the conduction electrons with low-energy scale below A^{x). In other words, 
the renormalization evolutions given by Eqs. (|6^)- (|65|) are stopped there. 



C. Scaling for Typical Cases 

Now we discuss whether or not the scaling properties, discussed in the preceding sections, 
are really measurable as a crossover phenomenon. 

The bare couplings and phonon energy are expressed in terms of the original parameters 
specifying our model as follows: 



(84) 



rs^ rs^ I I rsj 

D Ep \M 



b - I:, (85) 



where Ep is the Fermi energy of conduction electrons. Typical values of parameters are 
a ~ 10^^-10^^ and b ~ 10°. The crossover temperature Tk, (|79D , and the first-excited 
energy A^(Tk), (|82|) , of ion vibrations in the crossover region are expressed in terms of these 
parameters as 

— 1 3nab-\-l n , . 

Tk = 2~^a-^b^, (86) 
A^(Tk) = 2-(""^+i)a""^+^ (87) 



where we have assumed ^^(Tk) <^ 1 in Eq. (|82D. 

Tk and A^(Tk) depend sensitively on a and b, and their dependences are shown in 
Fig. 4 and 5, respectively. An anomalous behavior of physical quantities are observable for 
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sufficiently large Tk and small A^(Tk), which are realized for appropriate value oim/M and 
large one of g/D. It is because small value of m/M gives extremely small Tk, while large 
value of m/M gives large A^'(Tk) and small g/D provides small Tk and large A^(Tk). In 
other words, it is preferred that the electron-phonon coupling is large and the bandwidth is 
small. Such situation is expected to occur in d-band metals with narrow bandwidth such as 
A15 and C15 compounds p|,p9|. 

It is important to note that the region k-oT/D < A^{x) is not reached by the 
renormalization-group evolution if the parameters provide us with the crossover temperature 
Tk < A^{x). It is because the renormalization giving divergence of coupling constants f's 
are cut by the fictitious magnetic field acting on the pseudo-spin at around D'/D = A^{x). 

Typical cases of scaling are shown in Fig. 6. The case of 7k ^ A^(x) is shown in 
Fig. 6(a), where k^T/D, the dotted line, is always larger than A^(x), the solid line, in the 
right side of figure. In this case, anomalous behavior can be observed. It is noted that the 
coupling and increase rapidly below T < Ty^ while the coupling decreases reflecting 
the nature of two-channel Kondo effect. In the right side of figures (b) and (c), the energy 
scale quenching the pseudo-spin A determined by the intersection between the solid line, 
A^(x), and the dotted line, k-^T/D^ gives us the energy scale where the renormalization- 
group evolution should be stopped. The case of Tk < A is shown in Fig. 6(b). The 
ffist-excited energy of ion vibration is quenched before anomalous behavior sets in. The 
case of 7k ^ A is shown in Fig. 6(c). Any anomaly cannot be measurable because of 
extreamly small Tk. 

V. CONCLUSIONS AND DISCUSSIONS 

We have discussed the structure of renormalization-group evolutions of the strongly cou- 
pled electron-phonon system. The phonon degrees of freedom have been described by the 
Einstein oscillator. In order to investigate the low-energy physics, the Hilbert space of 
phonon has been restricted so that the states with more than two phonons are prohibited 



18 



though such states determine the algebra of the restricted phonon operators through the 
intermidiate or virtual states; and the conduction electrons are represented by spherical har- 
monics of k vectors on the Fermi surface because their relevant degrees of freedom is those 
around the Fermi level. 

The simplified Hamiltonian has the pseudo-spin degrees of freedom from the beginning 
due to phonons in the restricted Hilbert space, while the electrons have only a latent feature 
of the pseudo-spin which manifests itself through the renormalization steps. Namely, the 
scaling equations, which is derived by the multiplicative renormalization-group method at 
the two-loop level, shows that the polarization of conduction electrons corresponding to those 
of the phonon are selectively grown at the initial stage of renormalization-gruop evolution. 
Thus, the effctive Hamiltonian for low-energy physics are reduced to the anisotropic two- 
channel Kondo model with renormalized fictitious magnetic field acting on the pseudo-spin. 

In the anisotropic two-channel Kondo model with magnetic field, the scaling equations 
provide us with two characteristic energy scales, the crossover temperature Tk between the 
weak and strong coupling region and the renormalized first-excited energy of phonons 
which is the fictitious magnetic filed acting on the pseudo-spin. The anomaly associated 
with two-channel Kondo effect is observable when Tk > which is reahzed if the mass 
ratio M/m and the electron- phonon coupling g/ D are large enough. It is expected that such 
a sutiation is realized in the transition compounds with narrow bandwidth including A15 
and C15 compounds. 

The above result may be interpreted in terms of the adiabatic potential the ion feels as 
follows. As discussed in Refs. 2), 9), and 6), strong enough electron- phonon coupling makes 
the harmonic potential of individual Einstein oscillator be flattened and finally double-well 
shape in general. Namely, the adiabatic potential is deformed as shown schematically in 
Fig. 7. As a result, the first-excited energy of ion vibrations is softed considerably due to 
the two-channel Kondo correlation of pseudo spins. 

A salient feature of the above result is that the elcctron-phonon coupling suffers consid- 
erable renormalization due to pseudo-Kondo effect. Namely, the so-called Migdal approxi- 
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mation is not valid in this case because the higher order terms give logarithmically divergent 
contributions for the vertex correction. This latter situation may have been overlooked in 
the argument of proving a validity of the Migdal approximation. So the way that the Migdal 
approximation is broken down in the present theory appears to be somewhat different from 
that discussed by Varma [^], and Suzuki and Motizuki |^ in which the fact that the 
correlation length of CDW transition is very short is the origin of its breakdown. 

It is well recognized that the spin susceptibility is not renormalized by the electron- 



phonon interaction as far as the Migdal approximation is valid However, it seems 

still remain as a controversy whether the breakdown of the Migdal approximation directly 
implies the existence of renormalization of the spin susceptibility due to electron-phonon cou- 
pling P, p9| , p0| . The present result suggests that only the electron-phonon coupling cannot 
afford to renormalize the spin susceptibility which corresponds to the channel susceptibility 
in our fixed-point two-channel Kondo model. In order to obtain the enhancement of channel 
susceptiblity in that model, we need the perturbations breaking the particle-hole symme- 
try, such as the Coulomb repulsion among conduction electrons or potential scattering, as 
shown in Ref. 21) by the method of numerical renormalization group. In this sense, together 
with the Coulomb repulsion among conduction electrons, the interaction (P^, represent- 
ing the potential scattering, may be important to discuss a possible anomaly of the (real) 
spin susceptibility. It is intersting to note the logarithmic temperature dependence of spin 
susceptiblity observed in VsSi M is consistent with the scenario presented in this paper. 
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APPENDIX A: 



In this Appendix the solutions of equations (74)-(^) are given exactly and their simpli 
fied forms including only the sigular terms are derived for the case v^{xi) <^ v^{xi) <^ 1. 
First, Eq. ([75|) is rearranged to eliminate as 

dv^/dx 



.x\2 



4(1 - 2nv' 



Then substituting this into Eq. ([7^ ) multiplied by f^, we obtain 



1 d 


dv^ /dx 


d 




Adx 


(1 - 2nv^f 


dx 


1 - 2nv^ 



It is easily integrated to give the relation 



1 



dx = -- 



1 



dv' 



4 (1 - 2nv')[{v'y + 2nC^v' - C^^ 
where C is determined by the boundary condition at 

12 \ 1/2 



c 



[v'^lxiW — [v^(xi 



By integrating 



n 



1 — 2nv^{xi) J 
we obtain v^{x) as an implicit form: 

(1 - 2nv'y [v%xi)]^ + 2nCh'{xi) - 



(Al) 



(A2) 



(A3) 



(A4) 



2n^C^ + l , 
H : — , = In 



v')^ + 2nC^v' ~ [l-2nv'{xi)Y 

v' + nC^ - CVn^C^ + 1 v%xi) + nC^ + CVn^C^ + I 



SCVn^C^ + 1 



v 



(A5) 
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For the case v^{xi) ^ v^{xi) ~ v^{0) <^ 1, keeping only the singular terms in right-hand 
side of Eq. , we obtain 



X — Xi = — — ln[t>'*(xi)] — ^ In 



v^{xi) 



+ 



In 



v^{xi) 



+ 



In 



(A6) 



We can also obtain the explict form for xi from Eq. ( [72| ) as 



1 



Xi 



4f^(0)[l -nf^(O)] 



In 



2v^{xi) 



n 



4f^(0) 



In 



2v^{xi) 
f^(0) 



(A7) 



Then, eliminating a;i from Eqs. (|A6D and (|A7|) , we obtain Eq. (|77D for v^{x), and then we 
also obtain Eq. @ for v'^ix) from Eqs. (|AD and (0). 

Second, we eliminate dx from Eq. (^T]) by using Eq. (|A^). Then, the obtained equation 
can be integrated with the boundary condition at x = as follows: 



A^(0) 



1 — 2nv^ 
- 2nv^{0)^ 



1/2 



X 



X 



V' + nCl + Cjn^Cl + 1 i;^(0) + nCl - Cjn^Cl + 1 



nCo/2^n2c2+l 



(A8) 



where Cq is given by substituting v^{xi) and f^(xi) by f^(0) and f^(0), respectively, in 
Eq. ( [A4D . For the case f^'(O) <^ f^(0) ^ 1, keeping only the singular terms of the original 
couplings, we obtain the Eq. (|82|). 
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Figure Captions 

Fig. 1 

Diagrams for vertex correction of first order. Solid (dotted) line stands for conduction 
electron (psudo-spin) . 

Fig. 2 

Diagrams for vertex correction of second order. 
Fig. 3 

Diagram for Self-energy of pseudo-fermion of first order. 
Fig. 4 

{mjMfl^ dependence of Tk and A^(Tk). 
Fig. 5 

gjD dependence of 21< and A^(Tk). 
Fig. 6 

Typical cases of scaling; (a) anomalous behavior can be observed, (b) the first-excited 
energy of ion vibration is quenched before anomalous behavior sets in, (c) Tk is ex- 
tremely small. 

Fig. 7 

Schematic shape of adiabatic potential of ion; (a) for the case without electron-phonon 
coupling or the high-energy processes, and (b) for the fixed-point Hamiltonian or the 
low-energy-processes. 
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Fig. 1(a) 
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